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Abstract 

We have combined perturbative unitarity and renormalisation group equation arguments in order to 

m 

find a dynamical way to constrain the space of the gauge couplings (g[, g) of the so-called "Minimal 
we have shown that perturbative unitarity imposes an upper bound that is generally stronger than 



X 



Z' Models" . We have analysed the role of the gauge couplings evolution in the perturbative stability 
of the two-to-two body scattering amplitudes of the vector and scalar sectors of these models and 



the triviality constraint. We have also demonstrated how this method quantitatively refines the 
usual triviality bound in the case of benchmark scenarios such as the U(l) x , the U(1)r or the 
"pure" U(1)b-l extension of the Standard Model. Finally, a description of the underlying model 
structure in Feynman gauge is provided. 
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I. INTRODUCTION 



Nowadays the phenomenological importance of Beyond the Standard Model (BSM) physics 
at the TeV scale is recognised by the global experimental effort at the Large Hadron Collider 
(LHC). 

It is common belief that a Z' boson is among the first new objects that can potentially be 
detected at the LHC. The existing extensive literature is testimonial to the growing interest 
in them (see e.g. [l-4]). A particularly interesting class of theoretical scenarios incorporating 
a Z' boson are the so-called " (non-exotic / non- anomalous) Minimal Z' Models", extensively 



studied in the recent years 



io|. 



These models are based on an extension of the Standard Model (SM) gauge group with 
a further U(l) symmetry factor. The anomaly cancellation conditions imply the inclusion 
of three generations of right-handed neutrinos in the fermion sector, while the breaking of 
the new gauge group is provided by an extra singlet Higgs boson (thereby making the Z' 
boson a massive state). 

The purpose of this paper is to show that renormalisation group equation (RGE) based 



techniques 
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15| as well as a standard perturbative unitarity criterion 16( can be combined 



to give a dynamical way to constrain the two gauge couplings (g[ and g) of a set of Minimal 
Z' Models, with a particular attention devoted to some benchmark scenarios such as the 
"minimal" U(1)b-l, the U(1)r (no fermion charge associated to the left-handed fermions) 
and the S'O(10)-inspired U(l) x extensions (see |4| for an extensive overview). 

To this end, we propose a detailed study of the Goldstone and Higgs sectors of this kind of 
models with a view to extract the most stringent bounds on the (evolving) gauge couplings. 
We will make a comparison between this method and triviality arguments, showing that 
calling for perturbative unitarity stability conditions gives stronger constraints on g[ and g 
with respect to traditional triviality assumptions over most of the parameter space. For an 
exhaustive description of the theoretical setup and of our conventions see where also 
the RGE equations can be found. Finally, regarding perturbative unitarity techniques, we 
will expand below upon the methodology outlined in 18]. 

This work is organised as follows: in Section |TT] we introduce our parametrisation of 
the Scalar Lagrangian of the Minimal Z' Models, in Section II III we describe the theoretical 
methods adopted to constrain the gauge couplings, in Section [TV] we present our numerical 
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results while in the last section we give our conclusions; in Appendix [A] we discuss the 
gauge-fixing Lagrangian of Minimal Z' Models, in Appendix |B] we list the set of Feynman 
rules that is relevant in our calculation, and in Appendix O we give some explicit analytical 
results that have been used in this paper. 

II. THE PARAMETRISATION OF MINIMAL Z' MODELS 



We describe here our parametrisation of the Minimal Z' Models. Following [17|], the SM 
gauge group is augmented by a U(l) factor, related to the Baryon minus Lepton (B — L) 
gauged number. In the complete model, the classical gauge invariant Lagrangian, obeying 
the SU(3)c x SU{2)l x U{1)y x U(1)b-l gauge symmetry, can be decomposed as: 

= + + -^f + • (1) 

In this paper we are mainly interested in the scalar part of the Lagrangian 

% 8 = (D»H) ] D^H + (D» X ? D,X - V(H, X ) , (2) 
with the scalar potential given by 



V(H, X ) = m 2 H^H + f i 2 \x\ 2 + (H^H \ X \ 2 ) 



Ai f \ / H^H 



\2 J \\X 

= m 2 H^H + /i 2 | X T +Ai(# f #) 2 + A 2 | x I 4 +A 3 ^ f ^ | X ? , (3) 

where H and X are the complex scalar Higgs doublet and singlet fields, respectively. 

We generalise the SM discussion of spontaneous Electro-Weak Symmetry Breaking 
(EWSB) to the more complicated classical potential of equation ([3]). To determine the 
conditions for V(H, X ) to be bounded from below, it is sufficient to study its behaviour for 
large field values, controlled by the matrix in the first line of equation 03]). Requiring such 
a matrix to be positive-definite, we obtain the conditions: 

4AiA 2 - A 2 > , (4) 

A 1; A 2 >0. (5) 
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If the above conditions are satisfied, we can proceed to the minimisation of V as a function 
of constant Vacuum Expectation Values (VEVs) for the two Higgs fields. In the Feynman 
gauge, we can parametrise the scalar fields as 

V2 ^ v + (h + # 3 ) y V2 

where = 0i =F £02 are the would-be Goldstone bosons of W^, while 03 and 04 will mix 
to give z and z', the would-be Goldstone bosons of the Z and Z' bosons, respectively. The 
real and non-negative VEVs are v and x, for the Higgs doublet and singlet, respectively. 

We denote by hi and h 2 the scalar fields of definite masses, and m/i 2 respectively, and 
we conventionally choose m 2 ^ < mi ■ After standard manipulations, the explicit expressions 
for the scalar mass eigenvalues and eigenvectors are: 



X lV 2 + X 2 x 2 - y/{X lV 2 - X 2 x 2 ) 2 + (X 3 xv) 2 , (7) 



m 



h 2 



X lV 2 + X 2 x 2 + ^(X^ 2 - X 2 x 2 ) 2 + (X 3 xv) 2 , (8) 



cos a — sin a 
sin a cos a 




(9) 



where -f < a < f fulfils 1 : 



sin 2a = = , (10) 

a/ (Aif 2 - A 2 a; 2 ) 2 + (X 3 xvf 

o Aif 2 - A 2 x 2 

cos 2a = — . (11) 

^(X x v 2 - A 2 x 2 ) 2 + (X 3 xv) 2 

For our numerical study of the extended Higgs sector, it is useful to invert equations (J7|), 

(jHJ) and ffTU]) to extract the parameters in the Lagrangian in terms of the physical quantities 

m^, rrih 2 and a: 

2 2 

A x = !^(l-cos2a) + ^(l + cos2a), 
4ir 4?r 

2 2 

X 2 = ^(l-cos2a) + ^(l + cos2a), 

A 3 = sm2a . (12) 



1 In all generality, the whole interval < a < 2ir is halved because an orthogonal transformation is invariant 
under a — > a + it. We could re-halve the interval by noting that it is invariant also under a —> —a if we 
permit the eigenvalues inversion, but this is forbidden by our convention fn\ < . Thus a and —a are 
independent solutions. 
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In order to determine the covariant derivative, we must introduce JZym, in which the the 
non-Abelian field strengths therein are the same as in the SM whereas the Abelian ones 
can be written as follows: 



cv?Abcl 



'YM 



4 ^ 4 ^ 



(13) 



where 



UV 



daB v ~ d u B, 



(14) 
(15) 



In this field basis, the covariant derivative is: 



D^ = d u + ig s T a G« + + WiYB, + i(gY + g[Y B _ L )B' 



(16) 



To determine the boson spectrum, we have to expand the scalar kinetic terms like for 
the SM. As for the gauge bosons, we expect that there exists a mass-less gauge boson, the 
photon, whilst the other gauge bosons become massive. The extension we are studying is in 
the Abelian sector of the SM gauge group, so that the charged gauge bosons W ± will have 
masses given by their usual SM expressions, being related to the SU(2) l factor only. The 
gauge boson spectrum is then extracted from the kinetic terms in equation fl2]): 



= \d»hd a h + \{h + v) 2 (0 1) \gW/a a + g x B» + gB'» 

gauge Z o - 

= \d»hd a h + \{h + v) 2 \g 2 \W/ -iW/\ 2 

+ (gW/- gi B»-gB'»y 




(17) 



and 



gauge 



l -d»tid»ti + l -{ti + x) 2 {g[2B>») 2 , 



(18) 



where we have taken Y^ B L 
terms (see 



19 



2 in order to guarantee the gauge invariance of the Yukawa 
for details). In equation ( IT7|) we can recognise the SM charged gauge 
bosons W ± , with M\y = gv/2 as in the SM. The other gauge boson masses are not so 
simple to identify, because of mixing. In fact, in analogy with the SM, the fields of definite 
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mass are linear combinations of B^, W% and B ,fM . The explicit expressions are: 



( 



cos •&.„ — sin t?„, cos i?' sin sin i? 



V 



sin tf w cos i? w cos 
sin i?' 



cos sin 
cos i?' 



with -f < i?' < f , such that: 



and 



tan 2$' = 



^ + 16(f)V-^ 2 -^i 2 



V z *7 



(19) 



(20) 



M A = 0, 



M: 



where 



C = g 2 + gl + g 2 + iQ g? , 

D = 64v 2 x 2 (g 2 + g^g' 2 . 



(21) 



(22) 
(23) 



As for the Goldstone boson spectrum, it is possible to find a convenient way to write 
the mass matrix. Being H ~ (1,2,|,0) and x ~ (1? 1)0,2) the Higgs representations 
associated to each gauge group, in the gauge- Goldstone 2 bosons basis we find the following 
representation of the co-variant derivative: 



V = 



I 19 

o \g 










19 










\ 



(24) 



\ -\g -2xg[ J 

In the t'Hooft-Feynman gauge, it can be verified that the vector boson mass matrix is 
given by my = V(V) T . The related Goldstones mass matrix can as well be evaluated as 



rn 



(25) 



The 5x4 matrix follows from the five gauge bosons W l | j=i,3, Z, Z' and the four Goldstone bosons > 
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therefore we get 

/ ^g 2 \ 



rn 2 



o v 4g 2 o o 



•i ■ 

„2 



(26) 



^(^+^3+^2) XV gg< x 

\ ara&J 4x 2 (^)V 

The mass matrix in equation fl26|) shows that the Goldstones of the TV-boson have a 
mass that is equivalent to the SM one, while the 03 and 04 fields mix, as it happens for 
the related gauge bosons. We can diagonalise the neutral Goldstone block by means of a 
rotation of angle a g , defined by: 

tan 2a g = 5 , 27 

g 2 + g! + ^-16(^0 

obtaining, as expected, the neutral gauge boson masses as eigenvalues of the neutral Gold- 
stone boson sub-matrix. As for the neutral gauge boson sector, the Goldstones mix only if 
g 7^ 0. Finally, the neutral Goldstone bosons fulfil 



cos a g sin Otg 



sin a g cos a g 




(28) 



Now that the scalar Lagrangian has been presented in the Feynman gauge, we have all 
the elements to carry on with our analysis. Although not relevant for the latter, we also 
present for completeness the gauge-fixing Lagrangian in Appendix lAl 

The generic model that has been previously introduced spans over a continuous set of 
minimal U(l) extensions of the SM, that can be labelled by the properties of the charge 
assignments to the particle content. 

As any other parameter in the Lagrangian, g and g[ are running parameters, therefore 
their values have to be set at some scale. A discrete set of popular Z' models (see, e.g. 0, Q) 
can be recovered by a suitable definition of both g and g[. 

Even though we present results in the generic (g[-g) space, we will comment on a subset 
of particular interest: the "pure" B — L extension U(l) b-l ((Jew — 0) has a vanishing mixing 
between the massive neutral gauge bosons at tree-level at the EW scale, the 5'O(10)-inspired 
extension U{l) x (gEW = — 4/5^^) preserves the mixing ratio at any energy scale and the R 
minimal extension U(1)r (§ew — ~^g'i) in which a Z' is coupled to right-handed fermions 
only. In Table |T] we summarise these models emerging from our parametrisation. 
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Model 


Parametrisation 


U(1) B -L 


9EW = 


U(l)x 


9ew = -4/ 5 5i 


U(1)r 


9ew = -2g[ 



TABLE I: Specific parametrisations of the Minimal Z' Models: U(1)b-l, U(1) x and U(1)r. 
III. CONSTRAINING THE g[ - g SPACE 

Since it has been proven that perturbative unitarity violation at high energy occurs only 
in vector and Higgs boson elastic scatterings, our interest is focused on the corresponding 
sectors that have been already presented in Section HT1 



Following the Becchi-Rouet-Stora (BRS) invariance (see 211]), the amplitude for emis- 
sion or absorption of a "scalarly" polarised gauge boson becomes equal to the amplitude 
for emission or absorption of the related would-be-Goldstone boson, and, in the high en- 
ergy limit (s ^> vn? w± z z')i ^ ne amplitude involving the (physical) longitudinal polarisation 
(the dominant one at high energies) of gauge bosons approaches the (unphysical) scalar 



one, proving the so-called Equivalence Theorem (ET), see 22]. Therefore, the analysis of 
the perturbative unitarity of two-to-two particle scatterings in the gauge sector can be per- 
formed, in the high energy limit, by exploiting the Goldstone sector instead (further details 



of this formalism can be found in 18| ) 



Moreover, since we want to focus on g[ and g limits, we assume that the two Higgs bosons 
of the model have masses such that no significant contribution to the spherical partial wave 
amplitude (see below) will come from the scalar four-point and three-point functions (that 
is mi j2 <C 700 according to [isd), i.e. the Higgs masses are well below the Lee-Quigg-Tacker 



(LQT) limit 16] . It is important to remark that relatively high values of the Higgs masses, 
far below the unitarity limit, tend to lead to quartic coupling to values that become non- 
perturbative at high scales. On a side, this could considerably refine the unitarity bounds. 
On the other side, it could be non-consistent by triviality arguments (as a general rule, the 
larger the cut-off, the smaller the acceptable value of the Higgs mass). Beyond any doubt, 



given a cut-off energy, a good choice for the Higgs masses is the one explored in 17]. With 



this choice we exclude any other source of unitarity violation different from the size of the 



S 



g[ and g gauge couplings. 

Firstly, we focus on the techniques that we have used to obtain the aforementioned 
unitarity bounds in combination with the RGE analysis: for this, it is crucial to define the 
evolution of the gauge couplings via the RGEs and their boundary conditions. As already 



established in 



171 . |23| . the RGEs of g, g±, g[ and g are: 



d{g) 
d(logA) 

d(gi) 
d(logA) 

d(logA) 

d(g) 



167T 2 
1 

16vr 2 
1 

16vr 2 
1 



19 

"IP 

41 3 

6** 



41 



+ 2^ 2 ) + 2^W + rf) + 12<ffi 
6 3 



(29) 



d(logA) 16tt 2 

where g(.EW) ~ 0.65 and ^(.EW) — 0.36. This fully fixes the evolution of g[ and g with 
the energy scale A. 

In the search for the maximum g[(EW) and g(EW) values allowed by theoretical con- 
straints, the contour condition 



g[(A),g(A)<k, 



(30) 



also known as the triviality condition, is the assumption that enables one to solve the above 
system of equations and gives the traditional upper bound on the g[-g(EW) space at the 
EW scale. 

It is usually assumed either k = 1 or k = \/4~7T, calling for a coupling that preserves the 
perturbative convergence of the theory. Nevertheless, we stress again that this is an "ad 
hoc" assumption. Our aim, instead, is to extract the boundary conditions by perturbative 
unitarity arguments, showing that, under certain conditions, it represents a stronger con- 
straint on most of the gauge couplings parameter space. For this, we exploit the theoretical 
techniques that are related with the perturbative unitarity analysis, since they can be used 
to provide constraints on the theory, with a procedure that is not far from the one firstly 
described in detail by fl6|. 

A well known result is that, by evaluating the tree- level scattering amplitude of longitudi- 
nally polarised vector bosons, one finds that the latter grows with the energy of the process 
and, in order to preserve unitarity, it is necessary to include some other (model dependent) 
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interactions (for example, in the SM one needs to add the Higgs boson) and these must 
fulfil the unitarity criterion (a gain in the SM, the Higgs boson must have a mass bounded 
from above by the LQT limit 24J]). 

As already intimated, we also know that the ET allows one to compute the amplitude 
of any process with external longitudinal vector bosons Vl {V = W , Z, Z'), in the limit 
m v <C s, by substituting each one of these with the related Goldstone boson v = w ± , z, z' , 
and its general validity has been proven in {22 1. Schematically, if we consider a process with 
four longitudinal vector bosons, we have that M(VlVl — > VlVl) = M(vv — > vv) + 0(m v /s). 

While in the search for the Higgs boson mass bound it is widely accepted to assume small 
values for the gauge couplings and large Higgs boson masses, for our purpose we reverse such 
argument with the same logic: we assume that the Higgs boson masses are compatible with 
the unitarity limits and we study the two-to-two scattering amplitudes of the whole scalar 
sector, pushing the size of g[ and g up to the unitarisation limit. 

This limit is a consequence of the following argument: given a tree-level scattering ampli- 
tude between two spin-0 particles, M(s, 8), where 9 is the scattering (polar) angle, we know 
that the partial wave amplitude with angular momentum J is given by 

a J = ^T \ d(cos6)Pj(cosd)M(s,d), (31) 



where Pj are Legendre polynomials, and it has been proven (see 24J) that, in order to 



preserve unitarity, each partial wave must be bounded by the condition 

|He(aj(a))| < \ (32) 

By direct computation, it turns out that only J = (corresponding to the spherical 
partial wave contribution) leads to some bound, so we will not discuss the higher partial 
waves any further. 

Assuming that the Higgs boson masses do not play any role in the perturbative unitarity 
violation, we have verified that the only divergent contribution to the spherical amplitude 
is due to the size of the g[ and g couplings in the t-channel intermediate Z and Z' vector 
boson exchange contributions. In Appendix [B] we list the relevant 3-point Feynman rules 
that connect any of the two (external) scalars with either a Z or Z' (mediator). Hence, 
the relevant channels are represented by the 6-dimensional (symmetric) scattering matrix in 
Table HT1 plus the decoupled eigenchannel w + w~ — > w + w~ . 

10 





z z 


z z' 


z'z' 


h\hi 


hih 2 


/i2^2 


z z 

















z z' 

















z'z' 

















h\h\ 






r^i 











hih 2 

















hzh-i 


















TABLE II: Scattering matrix: we have used the simbol ~ just for illustrating the presence 
of a non-zero element in the correspondent scattering channels. 



After explicit evaluation, the spherical amplitude of the decoupled w + w eigenchannel, 
in the high energy limit, is: 



■J it 



16tt 



1 + 4 log 



Mz 
A 



+ 



£z' 

16tt 



1 + 4 log 



M z , 
A 



(33) 



and each element of the scattering matrix presents the following structure: 



S{Sj 



f*. 

167T 



1 + 4 log 



Mz 
A 



ho 



16n 



1 + 4 log 



M z , 
A 



(34) 



where S is a symmetry factor that becomes l/\/2 if the (initial or final) state has identi- 
cal particles, 1 otherwise, and A represents the scale of energy at which the scattering is 
consistent with perturbative unitarity, i.e. it is the evolution energy scale cut-off. It is im- 
portant to notice that the masses of the Z and Z' act as a natural regulator that preserves 
both the amplitude and the spherical partial wave from any t-channel collinear divergence 
and that both of them are completely defined by the choice of the gauge couplings and 
VEVs (see equation ( l21i) ). The non-vanishing coefficients of equations ( l33i) -( l34l) are listed 
in Appendix O 

It is well-known 3 that the most stringent unitarity bounds on the g[-g space are derived 
from the requirement that the magnitude of the largest eigenvalue of the scattering matrix 
does not exceed 1/2. 



3 The diagonalisation of the scattering matrix usually leads to stronger bounds not only in the S'M-case 
but also in BSM scenarios (e.g. |25j|V 
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Finally, if we consider the contour of this inequality, we find exactly the boundary con- 
ditions that solve the set of differential equations in ( 129]) . giving us the upper limits for g[ 
and g at the EW scale. In the next section we will combine all these elements to present a 
numerical analysis of the allowed domain of the gauge couplings. 



IV. RESULTS 



The set of differential equations (1291) has been integrated with the well-known Runge-Kutta 
algorithm and both the unitarity (equation (152"j) ) and triviality (equation f )5U]) ) conditions 
have been imposed as a two-point boundary value with a simple shooting method, that 
consisted in varying the initial gauge coupling values in dichotomous-converging steps until 
the bounds were fulfilled. 

Apart from the gauge couplings, other parameters play a role in the computation: the 
VEVs have been chosen in such a way that both Mz and Mz> are within the allowed 
experimental range (see 26( and 27J, respectively), and further that Mz> lies in the (9(1 — 10) 
TeV range, so that it does not spoil the high energy approximation Mz< "C A. By direct 
computation we verified that the Higgs mixing angle a does not play any significant role in 
the analysis, hence for each analysed point of the gauge couplings parameter space we have 
averaged the spherical wave greatest eigenvalue over the range — 1 < sin a < +1, finding 
a standard deviation never greater than ~ 2%. Finally, for illustrative purposes, we have 
chosen the triviality condition to be fixed by k = 1. 

As initial step of our numerical analysis, we have verified by direct computation that 
the spherical wave associated to the decoupled eigenchannel w + w~ — > w + w~ gives always a 
negligible contribution with respect to the greatest eigenvalue of the spherical wave scattering 
matrix in Table [III Therefore, in Figure [TJ we have overlapped the contour plots of both 
the greatest eigenvalue of the spherical wave scattering matrix allowed by unitarity and 
the gauge couplings allowed by triviality in the g[-g plane for the following values of the 
evolution/cut-off energy: A = 10 11 GeV (Figure Ha}, A = 10 15 GeV (Figure HE}, A = 10 19 
GeV (Figure [Tc]). It is clear that the boundary condition imposed by the perturbative 
unitarity stability (dashed lines) constrains the parameter space considerably more than the 
well-known triviality bound (dotted lines). In few cases the triviality bound is (slightly) 
more important than the unitarity bound: this condition is realised at energies <C 10 19 
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(c) 
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(d) 



FIG. 1: Contour plot of both the greater eigenvalue of the spherical wave scattering matrix 
(Table [TTJ) allowed by unitarity (dashed lines) and the gauge couplings allowed by triviality 
(dotted lines) for several values of the cut-off energy: A = 10 11 GeV (light green/grey lines, 

Figure [Taj, A = 10 15 GeV (dark red/grey lines, Figure Hbj) . A = 10 19 GeV (black lines, 
Figure [Tc]). Figure [Tdl shows a summary of the most stringent bounds at different values of 
the cut-off energy, with focus on some peculiar parametrisation of g (Table H]). 
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Lo 5l0 (A/GeV) 


1 


9 


11 


1 o 

lo 


15 


1 ^7 
1 / 


1 n 
19 




l 


0.594 


0.527 


0.477 


0.439 


0.407 


0.380 


0.357 


r/Yi \ 
u ( l )b-l 




















U 


n A Q7 


U.oyU 


U.OOO 


u.zy / 


U.zoy 


n o a v 
U.Z4W 


U.zzy 




T 


0.682 


0.620 


0.573 


0.535 


0.504 


0.479 


0.457 


U(l)x 




















U 


0.531 


0.424 


0.364 


0.324 


0.295 


0.272 


0.254 




T 


0.362 


0.328 


0.300 


0.276 


0.254 


0.235 


0.218 


U(1)r 




















U 


0.429 


0.344 


0.293 


0.258 


0.232 


0.210 


0.192 



TABLE III: Triviality bounds (with k = 1) and unitarity bounds on g[ in (non-exotic) 
Minimal Z' Models, for several values of the energy scale A. 

GeV and g = hg[ where \h\ > 2 (see Figure llatil b|) . Otherwise, the unitarity condition 
noticeably refines the bounds on the allowed parameter space considerably, as it is clear 
from Figure Hd] in which the most stringent bounds are plotted for the aforementioned 
values of the cut-off energy. In the same figure, we plotted three lines as reference for some 
peculiar parametrisation of g already mentioned in Section [III it is clear that for each one 
of these models the unitarity condition is always more important than the triviality one. 

As for these specific parametrisations, in Figure [2] we have plotted the boundary value of 
g[ against the evolution/cut-off scale A, using both the perturbative unitarity stability con- 
dition (dashed lines) and the triviality condition (dotted lines). For the U{1)b-l (Figure l2a|) 
and the U{l) x (Figure l2b|) extensions of the SM model, the unitarity bound is always more 
stringent than the triviality one. For the U(1)r (Figure [2c]) extension, this is only true if 
A > 10 10 GeV. In Figure [2d] we have plotted the best bound on g[ (and then g) against the 
evolution/cut-off energy scale A. 

In order to summarise these results, in Table IIHI we present a comparison between the 
triviality and the unitarity bounds on g[ for several values of the energy scale A for our 
choice of Minimal Z' Models. 

V. CONCLUSIONS 

In this paper, we have shown that, by combining perturbative unitarity and RGE methods, 
one can significantly constrain the gauge couplings (g[ and g) of a generic/universal (non- 
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FIG. 2: The g[ (g) couplings bounded either by triviality (dotted lines) or unitarity 
(dashed lines) conditions plotted against the evolution/cut-off energy for several peculiar 
choices of the gauge couplings parametrisation: U{1)b-l (black lines, g = 0: Figure l2al) . 
[/(l);^ (dark red/grey lines, g = —4/5g[: Figure l2bl . U(1)r (light green/grey lines, 
g — —2g[: Figure l2c|) . Figure l2dl shows a summary of the most stringent bounds at 
different values of the considered parametrisations of g (Table H]). 
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exotic/non-anomalous) Z' gauge boson, by imposing limits on their upper values that are 
more stringent than standard triviality bounds. (Also notice that, as unitarity is more 
constraining than triviality, the stability of the perturbative solutions obtained through the 
former is already guaranteed by the latter.) 



The present work, alongside 



Q, 



171 ] and (28j, is in particular part of the long-term 
effort to establish the theoretical bounds on the parameter space of the B — L based U(l) 
extension of the SM. 
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Appendix A: Gauge-fixing Lagrangian of Minimal Z' Models 



As for the Goldstone bosons sector, the mass matrix and interactions for the ghost fields are 
defined by the matrix V, as in equation (I24p . via 



m 



ghost 



V{V) T . 



(Al) 



Notice that the m 2 ghost and the ml of equation ( 1251) have different numbers of zero- 
eigenvalues, but their non-zero eigenvalues are in a one-to-one correspondence; furthermore, 
the eigenvalues of the gauge- fixing mass matrix are the same of the gauge boson mass matrix. 

Then, the ghost Lagrangian is defined, in the t'Hooft-Feynman gauge, as 



ghost 



-c a 



{d^D^ + V a ■ (V b + S b ) 



(A2) 



where the matrix S represents the link between the fluctuations (Goldstones) of the Higgses 
around their VEVs and the gauge bosons; a convenient way to write this matrix is 

2 

9 



v 



9 -h 

2 


9 A 


9 A 


91 A 


2 h 




9 A 


-h 
2 


2 n 


2 Y 




9 A 




9 -h 
2 


2 




T 















>1 



9 



h 



(A3) 



-2h'g[ J 



Finally, the ghost fields ( c ) read as 

w{ w 9 2 w 9 3 B9 (B'Y 



(A4) 
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Appendix B: Feynman rules associated with a neutral gauge boson exchange in a 
scalar two-body scattering 



In the following we list the set of Feynman rules that enter in the calculation described 
Section 1111} these have been obtained by means of implementing the information of the 
scalar Lagrangian (see Section HT1 and Appendix |A"1) in the LanHEP package [29]; all the 
momenta p's are considered in-coming: 

h x - Z-z :=> — (s w s g > c w c a c ag g p£- s w s e , c w c a c ag g p^ 



+ c a c ag Cqi e p^ - c a c ag c i eVh~ ^ s w s a s ag s e > c w g 1 pi 
+ As w s a s ag s e i c w g[ p%) (Bl) 
hi- Z - z' :=> - — (s w s ag s e , c w c a gp^- s w s ag s e , c w c a g p%, 

ZC\v Sw 

+ s ag c a cqi e p%, - s ag c a cq> e p£ + As w s a s > c w c ag g[ p^ 

- As w s a sqi c w c ag g[ p^ z ) (B2) 
hi - z - Z' :=> - — (s w c w c a c ag cg< g pfc- s w c w c a c ag c e > g p^ 

lCw Sw 

- se' c a c ag + sg> c a c ag epjj- 4s w s a s ag c w c 6 > g[ p^ 

+ As W S a S ag C W Cqi g[ pt[) (B3) 
hi — Z' — Z' (s W Sag C W C a Cqi g p^ - S W S ag C W C a Cqi g <p%, 

Zcw sw 

- Sag Sgi C a C p£ + S ag Sg> C a C p^ h + 4s W S a C W C ag Cqi g[ p^ 

- As W S a C W C ag Cqi g[ p£) (B4) 
h 2 - Z - Z :=>• (s W S a Sqi C W C ag fir Ph - S W S a Sqi C W C ag fir 

Zcw Sw 

+ S a Cag Cqi C p^ - S a C ag Cqi C p^ h + As W S ag Sqi C W C a g[ p£ 

- As w Sag sqi c w c a g[ p£) (B5) 

hi- Z-z' :=> (s W S a Sag Sqi C W fi Ph ~ S W S a S ag Sqi C W g p£ 

lCw Sw 

+ S a Sag Cqi C p£, - S a S ag Cqi Cp^- 4s W Sqi C W C a C ag g[ P% 

+ As w sqi c w c a c ag g[ P z i) (B6) 
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hv-z-Z' — (s w s a c w c ag c e , g pfc- s w s a c w c ag c e , g p^ 

- s a s e > c ag ep^ + s a s e > c ag ep{J + As w s ag c w c a cg> g[ p^ 

- 4s w s ag c w c a c e > g[ p£) (B7) 
h 2 - Z' - z :=>• - — (s w s a s ag c w eg, g p£ - s w s a s ag c w eg, g p^ 

ZCw Sw 

- s a s ag s e > e p^ + s a s ag s e > ep^- 4s w c w c a c ag cg> g[ p^ 

+ As w c w c a c ag c e > g[ p%) (B8) 



'I 

w + -w~ - Z :=> ((1 - 2s w ) c e , e + % s > c w g p^_ 

ZCw Sw 

- (1 - 2s w 2 ) eg, e p^ + - s w s e > c w g P^+) (B9) 

w + - w~ - Z :=>• — ((1 - 2s w ) s e > e - s w c w cg> g p^_ 

lew s w 

- (1 - 2s w ) s e , e p" w+ + s w c w cgi g p^ + ) (BIO) 
In the previous formulae, we have used the following notation: 

cw(sw) — > cos9 w (sm9w), 

c a (s a ) cos a(sin a), 
c a g{s ag ) ->■ cosa 9 (sina s ), (Bll) 

cg'(sgi) — > cos 9' (sin 9'), 

_^ 99i 
e — > 



\/9 2 + 9 2 i 
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Appendix C: Explicit value of the J*'? coefficients of equations (I33j) - (|34|) 

The non- vanishing coefficients related to the structure of equation (1341) . in the high energy 
limit, for each entry of the scattering matrix are the following: 

f z = 

J zz,h\h\ 

= \ {- c l c lg c e> {9 2 + 9 2 i) + {c a c ag g - 4:g'is a s a g) 2 sl) , (CI) 

J zz,hxhi 

= J (lQc 2 e ,(g' 1 )' 2 's 2 a s 2 ag - %c o ,c 0l gCQ l g' 1 s o ,s ag (c e ,g+ \J ' g 2 + gfsg^j 

+ c 2 a c 2 Qg (clg 2 + 2c e ,\Jg 2 + g 2 gs e , + (g 2 + gf) 4)) , (C2) 



f z 

J zz,h\ 



hi 



- ( Acagg^aSagSg, [ cg> \J g 2 + g{ - gs e > 



+ ^Cagg^agSo' I -cg> \J g 2 + g{ + gs e > 

+ c a s a f-16(^) 2 4 ff 4 + 4 9 (<$> (g 2 + gf) - ^\Jg 2 + glgs e > +g*4 jj J , (C3) 
f z ' 

= \ \ 4c 2 a c ag co>g[s ag (^ce>g+ \J ' g 2 + ~gfs e A 



- ^c a gCe,g x s a s a g \c e <g + \J g 2 + g{s e , 
+ c a s a ^-16c 2 ,(g[) 2 s 2 ag + c 2 ag (^c 2 ei g 2 + 2c e > \J ' g 2 + ~$gs v + (g 2 + gf) 4^^) > 



(C4) 



f z 

J zz,hih2 
1 



clgcl {g 2 + gl) 4 - 2c 2 ag c el \J g 2 + g\gs 2 a s e , 
- 8c a c ag c e/ \Jg 2 + g 2 g[s a s a gSe' + (c ag gs a + 4c 



(C5) 



f z = 

J zz,h,2h,2 

4 {cg>(c ag gs a + 4:C a g[s ag ) 2 - c 2 ag (g 2 + gf) s 2 a s 2 gi ) , (C6) 
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ez 

J zz',hihi 
1 



±6c ag (g[) 2 s 2 a s ag s 2 9/ + 4c a g[s a (c 2 ag + s 2 ag ) s e , (c e >\j 9 2 + g\ - gse^j 
+ c 2 a c ag s ag (c 2 e , (g 2 + gl) - 2c e ,^g 2 + gfgs e > + <? 2 4^) ^ , (C7) 



fzz' 



zz' ,h\hi 



i ^16c a9 c^(^) 2 s 2 Sa9 _ ^ Ca c el g' 1 s a (c 2 ag + s 2 ag ) (c e ,g + \Jg^ 2 ~Vg\s f)l 

[clg 2 + 2c e , y/g* + ^ + (<? 2 + 4) 4^ , (C8) 



+ c a c ag s ag I c 



fzz' 



zz',h\h2 



-co' \lg 2 + g\ + gs e > 



= \ (^g'islslgSe, (^ce'\J g 2 + gl~ gse^J + ^c 2 a c 2 ag g' lSe/ ^- 
+ c a c ag s a s ag ^4 (# 2 + # 2 ) - 2c e > \jg 2 + gfgs e > + (-16(g[) 2 + g 2 ) 4^ , (C9) 



fzz' 



zz' ,h\h2 



= \ (^c 2 a c 2 ag c e >g[ (^ce'9 + \j ' g 2 + gfse^J - ^ g' lS 2 a s 2 ag (^c 9 'g + \J ' g 2 + gjso^J 

+ c a c ag s a s ag (c 2 e , (-16(^) 2 + £ 2 ) + 2c e ,^g 2 + g 2 gs , + (4 + g\) 4^ , (CIO) 



ce>\l g 2 + g\- gs e > 



= J ^4c a ^s a 4 9 s e' (^ce' \fg 2 +~gl - gs e ^j - 4c a c 2 9 ^s a s e / 
+ c Qff s aff (4 (4 + 4) 4 - + 4s4^' + (-16c 2 (^) 2 + ^4) 4) ) , (CH) 



f z . 



zz',h,2ll2 



= \ (-16c 2 a c ag c 2 9/ (g[) 2 s ag + AcaCvg'^^clg - s 2 ag ) (c e ,g + \J g 2 + gls e ^j 

+ c ag s 2 a s ag (c 2 e ,g 2 + 2c e ^ g 2 + glgs 6 ' + (g 2 + g\) ^ , (C12) 
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fZ 

J z'z' ,h\h\ 
1 

4 



44 {g 2 + gl) s 2 ag + 8c a c ag c 6/ \J g 2 + g 2 g[s a s ag sg> 
- 2c lce'\Jg 2 + gfgs 2 ag s 9 ' + (-Ac ag g[s a + c a gs ag f s 2 e }j , (C13) 
f z ' 

J z'z' ,hih± 

= ^ (^§C 2 a gC 2 e'(g'l) 2 S 2 a ~ SCaCagCe'g'iSaSag (cg,g + \J ' g 2 + g^Sg^j 

+ c 2 a s 2 ag [clg 2 + 2c e , sjg 2 + g 2 gs e , + (g 2 + g\) s 2 ^, (C14) 

= fz'z' 



z'z' ,h\h,2 



+ ^CaCagg^ugSg, I -Cg/ \J g 2 + g{ + gSg> 

+ c a s a (4 (g 2 + gf) s 2 ag - 2c e/x Jg 2 + g 2 gs 2 ag se' - (l6c^(^) 2 - tfs 2 ^) s 2 }j ^ , (C15) 
f z ' 

J z'z',hih,2 

= J (^clCagCgig'^ag (cg,g + \J g 2 + g^g^j 



- Ac ag cg l g' 1 s 2 o ,s 0ig \ c e ,g + \J g 2 + g\sg> 

+ c a s a (-IScl/M) 2 + 4 (4,g 2 + 2co,yJg 2 + gfgs e , + (g 2 + gf) 4jjj , (C16) 

P = 

J z'z' \h2h2 

= \ ( c l {a 2 + gf) s l s l g - 2c 9' \/g~ 2 +~gls a s ag (4:C a c ag g' 1 + gs a s ag )sg> 

+ (^CaCa^ + gSaSag^Sg,), (C17) 

p' 

J z' z' ,h 2 h2 

= ~ U644 g 4(fl'i) 2 + 8c a c a9 c e /5fJs Q s a;? ^c e /<? + \fg 2 +~g 2 ~sg^j 

+ s 2 a s 2 a g (4s 2 + 2c e , ^ 2 + + (s 2 + gl) 4 J J • (C18) 

The non- vanishing coefficients related to the structure of equation ( 1331) . in the high energy 
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limit, are the following: 



f z = 

4 (9 2 - 9lf VF+gl + 2c e , (fl 4 - gf) gs e , + (g 2 + gff 2 ggj, 

4 {g 2 + gif 2 

fi 



*. (9 2 + 9if 2 Q 2 - 2c e , (g* - gf) gs e . + (g 2 - glf 



4(g 2 + g!f 2 

In the previous equations, we have used the notation of equations (IB lip . 



23 



[1] T. G. Rizzo, (1996), arXiv:hep-ph/9612440. 

[2] M. Cvetic and P. Langacker, (1997), arXiv:hep-ph/9707451. 

[3] A. Leike, Phys. Rept. 317, 143 (1999), arXiv:hep-ph/9805494. 

[4] M. S. Carena, A. Daleo, B. A. Dobrescu, and T. M. P. Tait, Phys. Rev. D70, 093009 (2004), 
arXiv:hep-ph/0408098. 

[5] T. Appelquist, B. A. Dobrescu, and A. R. Hopper, Phys. Rev. D68, 035012 (2003), arXiv:hep- 
ph/0212073. 

[6] P. H. Chankowski, S. Pokorski, and J. Wagner, Eur. Phys. J. C47, 187 (2006), arXiv:hep- 
ph/0601097. 

[7] A. Ferroglia, A. Lorca, and J. J. van der Bij, Annalen Phys. 16, 563 (2007), arXiv:hep- 
ph/0611174. 

[8] P. Langacker, Rev. Mod. Phys. 81, 1199 (2009), arXiv:hep-ph/0801.1345. 
[9] J. Erler, P. Langacker, S. Munir, and E. R. Pena, JHEP 08, 017 (2009), arXiv:hep- 
ph/0906.2435. 

[10] E. Salvioni, G. Villadoro, and F. Zwirner, JHEP 11, 068 (2009), arXiv:hep-ph/0909.1320. 
[11] A. D. Linde, JETP Lett. 23, 64 (1976). 
[12] A. D. Linde, Phys. Lett. B62, 435 (1976). 
[13] S. Weinberg, Phys. Rev. Lett. 36, 294 (1976). 
[14] K. G. Wilson, Phys. Rev. B4, 3184 (1971). 
[15] K. G. Wilson and J. B. Kogut, Phys. Rept. 12, 75 (1974). 
[16] B. W. Lee, C. Quigg, and H. B. Thacker, Phys. Rev. D16, 1519 (1977). 
[17] L. Basso, S. Moretti, and G. M. Pruna, Phys. Rev. D82, 055018 (2010), arXiv:hep- 
ph/1004.3039. 

[18] L. Basso, A. Belyaev, S. Moretti, and G. M. Pruna, Phys. Rev. D81, 095018 (2010), arXiv:hep- 

ph/1002.1939. 
[19] E. E. Jenkins, Phys. Lett. B192, 219 (1987). 

[20] W. Buchmuller, C. Greub, and P. Minkowski, Phys. Lett. B267, 395 (1991). 
[21] C. Becchi, A. Rouet, and R. Stora, Annals Phys. 98, 287 (1976). 
[22] M. S. Chanowitz and M. K. Gaillard, Nucl. Phys. B261, 379 (1985). 



24 



[23] L. Basso, A minimal extension of the Standard Model with B — L gauge symmetry, Master's 
thesis, Universita degli Studi di Padova, 2007, http://www.hep.phys.soton.ac.uk/~l.basso/B- 
L_Master_Thesis.pdf. 

[24] M. Luscher and P. Weisz, Nucl. Phys. B300, 325 (1988). 

[25] S. Kanemura, T. Kubota, and E. Takasugi, Phys. Lett. B313, 155 (1993), arXiv:hep- 
ph/9303263. 

[26] Particle Data Group, K. Nakamura et al, J. Phys. G37, 075021 (2010). 

[27] G. Cacciapaglia, C. Csaki, G. Marandella, and A. Strumia, Phys. Rev. D74, 033011 (2006), 

arXiv:hep-ph/0604111. 
[28] L. Basso, S. Moretti, and G. M. Pruna, (2010), arXiv:hep-ph/1009.4164. 
[29] A. V. Semenov, (1996), arXiv:hep-ph/9608488. 



25 



